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ABSTRACT 



Aims. We derive the evolution of the infrared luminosity function (LF) over the last 4/5ths of cosmic time, using deep 24 and 70 pm 
I ■ imaging of the GOODS North and South fields. 

Methods. We use an extraction technique based on prior source positions at shorter wavelengths to build the 24 and 70 /jm source 
catalogs. The majority (93%) of the sources have a spectroscopic (39%) or a photometric redshift (54%) and, in our redshift range of 
interest (i.e., 1.3 < z < 2.3) ~ 20% of the sources have a spectroscopic redshifts. To extend our study to lower 70 yum luminosities we 
perform a stacking analysis and we characterize the observed L 2 4/(i +z) vs L 7 o/(i +; ) correlation. Using spectral energy distribution (SED) 
templates which best fit this correlation, we derive the infrared luminosity of individual sources from their 24 and 70 pm luminosities. 
We then compute the infrared LF at z ~ 1.55 ± 0.25 and z ~ 2.05 ± 0.25. 

Results. We observe the break in the infrared LF up to z ~ 2.3. The redshift evolution of the infrared LF from z = 1.3 to z = 2.3 is 
consistent with a luminosity evolution proportional to (1 + z) L0±0 ' 9 combined with a density evolution proportional to (1 + z)~ L1±L5 . 
At z ~ 2, luminous infrared galaxies (LIRGs: 10 u L o <L IR < 10 12 L o ) are still the main contributors to the total comoving infrared 
luminosity density of the Universe. At z ~ 2, LIRGs and ultra-luminous infrared galaxies (ULIRGs: 10 12 L o < Lir) account for ~ 49% 
and ~ 17% respectively of the total comoving infrared luminosity density of the Universe. Combined with previous results using the 
same strategy for galaxies at z < 1.3 and assuming a constant conversion between the infrared luminosity and star-formation rate 
(SFR) of a galaxy, we study the evolution of the SFR density of the Universe from z = to z = 2.3. We find that the SFR density of 
the Universe strongly increased with redshift from z = to z = 1.3, but is nearly constant at higher redshift out to z = 2.3. As part of 
the online material accompanying this article, we present source catalogs at 24 ^m and 70 yum for both the GOODS-North and -South 
fields. 
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C3 ; 1- Introduction (SED) libraries (Chary & Elbaz 2001; Lagache et al. 2003; 

Dale & Helou 2002). All of these studies concluded that the 

The important contribution of infrared luminous galaxies relative contribution of ULIRGs to the SFR density of the 

(Luminous Infrared Galaxies, LIRGs: 10 n L o <Lir< 10 12 L o ; Universe increases with redshift, and may even be the dominant 

Ultra-Luminous Infrared Galaxies, ULIRGs 10 12 L Q <L IR ) in component at z ~ 2. However, these conclusions still need to 

the evolution of the star-formation rate (SFR) history of the be confirmed since there are large uncertainties at high-redshift 

Universe is now well established up to z ~ 1 (Chary & Elbaz in transforming observed 24 pm flux densities to far-infrared 

2001; Franceschini et al. 2001; Xu et al. 2001; Elbaz et al. luminosities (Papovich et al. 2007; Daddi et al. 2007a). To 

2002; Metcalfe et al. 2003; Lagache et al. 2004; Le Floc'h study the high-redshift evolution of the SFR density, one has 

et al. 2005; Magnelli et al. 2009). Their contribution to the to combine deep mid- and far-infrared observations in order 

SFR density of the Universe increases with redshift up to z ~ 1 to infer robust bolometric corrections and to clearly constrain 

where the bulk of the SFR density occurs in LIRGs. Study of the location of the break of the infrared luminosity function (LF). 
this evolution was made possible through the use of large and 

accurate spectroscopic and/or photometric redshift catalogs as At z ~ 2, the observed 70 pm emission corresponds 

well as deep 24 and 70 yum surveys obtained by Spitzer. approximately to the rest-frame 24 pm luminosity, which 

At z > 1.3, the SFR history of the Universe has been was proven to be a good SFR estimator in the local Universe 

derived by several studies (Perez-Gonzalez et al. 2005; Caputi (Calzetti et al. 2007). The reliability of this SFR estimator 

et al. 2007) using deep 24 pm imaging and infrared bolometric seems to hold at high-redshift since the SFR of z ~ 2 galaxies 

correction estimated from local spectral energy distribution estimated from their observed 70 pm flux densities and their 
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Fig. 1. (Left) Comparison between the photometric and spectroscopic redshifts of our 24 /jm-selected catalog. This comparison is made using 1670 
galaxies which have both kinds of redshifts. Dashed lines represent the relative errors found in the redshift range of our study (i.e., o" A -/ (1+ .) = 0.14 
for 1.3 < z <2.3). (Right) Uncertainty in determining L IR from 24 pm fluxes densities and the Chary & Elbaz (2001) library due to error in 
photometric redshift estimates (black histogram). The dashed line is a Gaussian fit to this uncertainty distribution with a = 0.15. Note that the 
uncertainty distribution is not Gaussian. 



radio emissions are in good agreement (Daddi et al. 2007b). 
Thus, to get robust estimates of the SFR of distant galaxies, we 
decided to use deep 70 pm observations obtained by Spitzer. 

The main difficulty of using 70 pm observations to study 
star-formation at z ~ 2 resides in the limited depth of the 
existing Spitzer data, even from the very deepest observations 
such as those in the Great Observatories Origins Deep Survey 
(GOODS). In this study, we overcome this limitation by using 
a stacking analysis. As shown in Papovich et al. (2007) this 
analysis allows characterization of the 24 vs 70 pm correlation 
and thus constrains the bolometric corrections to be applied 
to the 24 pm flux densities. Using deep 24 and 70 pm images 
of the GOODS-North and South fields we find that the 24 vs 
70 pm correlation observed at high redshift is significantly 
different from predictions by standard SED libraries. This 
deviation can be interpreted as a possible signature of an 
obscured active galactic nuclei (AGN) (Daddi et al. 2007a) or 
simply as a SED evolution characterized by stronger polycyclic 
aromatic hydrocarbon (PAH) emission (Papovich et al. 2007). 
Both interpretations are discussed and two different bolometric 
corrections are inferred. Based on these bolometric corrections 
we derive the infrared LF in two redshift bins (i.e., 1.3 < z < 1.8 
and 1.8 < z < 2.3). For the first time these infrared LFs take 
into account the evolution of SED observed in high-redshift 
galaxies. By comparison, extrapolation from the observed 
24 pm emission alone results in significant overestimates in 
their infrared luminosity (Papovich et al. 2007; Daddi et al. 
2007a) while extrapolations from the observed 850 pm emission 
constrains only the most luminous ULIRGs. 

Throughout this paper we will use a cosmology with 
H = 70 km s _I Mpc\Q A = 0.7 and Q M = 0.3. 

2. Data 

2. 1 . Infrared Imaging 

The 24 pm imaging of the GOODS-N (\2 h 36 m , +62°14') and 
GOODS-S (3 /! 32 ffl , -27°48') fields were obtained as part of 
the GOODS Legacy program (PI: M. Dickinson). The 70 pm 
data in both GOODS fields were obtained by Spitzer GO 
programs GO-3325 and GO-20147 (PI: Frayer). In the north 
they cover a region of roughly 10' x 16', while in the south 
they cover a somewhat smaller fraction of the GOODS 24 pm 



area, roughly 10' X 10'. The Frayer data have been combined 
with additional 70 pm observations covering a wider area 
from the Far-Infrared Deep Extragalactic Legacy (FIDEL) 
program (PI: Dickinson), as well as shallower data from Spitzer 
Guaranteed Time Observer (GTO) programs (PI: Rieke). While 
our catalog is complete to the observational limits, here we 
restrict our analysis to the deepest regions at 70 pm covered by 
the Frayer+FIDEL+GTO data, covering areas of 194 arcmin 2 
and 89 arcmin 2 in GOODS-N and -S, respectively. 

At the resolution of Spitzer at 24 and 70 pm, all sources in 
our fields are point sources (i.e FWHM~ 5.9" and 18"at 24 pm 
and 70 pm respectively). Flux densities at these wavelengths 
are hence estimated using a PSF fitting technique based on the 
knowledge of the expected positions of the sources (Magnelli 
et al. 2009). For the 24 pm data, we use the position of the 
IRAC 3.6 pm sources as priors. This choice is motivated by the 
fact that the IRAC 3.6 pm data are 30 times deeper than our 
current 24 pm observations and that the typical S 24/jm/S3.6/jm 
ratio spans the range [2-20] (Chary et al. 2004). Hence we can 
assume that all 24 pm sources have an IRAC 3.6 pm counterpart. 
For the 70 pm data we use as prior the IRAC positions of our 
flux limited sample of 24 pm sources. At this wavelength our 
choice is straightforward since the typical S jo/jm/S i\ m ratio 
spans the range [2 - 100] (Papovich et al. 2007) and the 24 
pm observations are about 100 times deeper than our current 
70 pm observations. These assumptions are tested by visual 
inspection of residual images which would reveal any 24 pm or 
70 pm sources missed due to the lack of priors. We find no such 
sources. 

For GOODS-N, we use the IRAC catalog generated from 
the publicly available GOODS Legacy data (19437 objects 
detected at 3.6 pm with a 50% completeness limit of 0.5yuJy). 
The GOODS-S IRAC data have been incorporated into the 
SIMPLE IRAC Legacy Survey observations covering the wider 
Extended Chandra Deep Field South, also observed with MIPS 
as part of the FIDEL Legacy program. We use the SIMPLE 
IRAC catalogs (Damen et al. 2011, 61233 objects detected at 
3.6 pm with a 50% completeness limit of 1.5 pSy) as priors 
for the MIPS source extraction over the whole FIDEL area, 
although here we concentrate only on the deep 10' x 10' 
GOODS-S region. 

Using Monte carlo simulations we are able to estimate the 
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quality of our 24 and 70 pm catalogs (see Magnelli et al. 2009). 
We find that, in both fields, the 24 and 70 pm observations reach 
an 80% completeness limit of 30 pJy and 2.5 mJy respectively. 
In the GOODS-N and S fields we detect 2151 and 870 sources 
respectively with S24 > 30 //Jy, and 119 and 50 sources with 
S 70 > 2.5 mJy. Tables B.8, B.9, B.10 and B.ll of the online 
material give excerpt of the complete GOODS-N/S 24 pm and 
70 pm catalogs that are now available on CDS 1 . These catalogs 
extend below the 80% completeness limit, and covers the full 
extent (approximately 10' X 16') of the GOODS-S region, 
not only the smaller 10' x 10' region with the deepest 70 pm 
imaging that is used for the analysis in this paper. 

Since calibration factors taken to generate the final 24 
and 70 pm mosaics are derived from stars, color corrections 
(at most ~ 10%) have to be apply to all our fluxes. These 
color corrections being highly dependent on the redshift and the 
intrinsic SED of the sources, we decided to introduce these color 
corrections directly in the ^-correction used to estimate the LF 
since both quantities are taken into account in this computation. 
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2.2. Redshifts 

In this study we use spectroscopic redshifts coming from a 
combination of various studies (Cohen et al. 2000; Wirth et al. 
2004; Cowie et al. 2004; Le Fevre et al. 2004; Mignoli et al. 
2005; Vanzella et al. 2006; Reddy et al. 2006; Barger et al. 2008; 
?, Kurk et al. in prep for GMASS redshifts and finally Stern et 
al. in prep). Photometric redshifts are computed using Z-PEG 
(Le Borgne & Rocca-Volmerange 2002) and all optical and near 
infrared data currently publicly available. In GOODS-N, optical 
observations in the BVIz passbands were obtained with the 
Advanced Camera for Surveys (ACS) onboard the Hubble Space 
telescope (HST) as part of the GOODS ACS Treasury program 
(M. Giavalisco and the GOODS Team, 2010, in preparation) 
while near infrared observations in the JK passbands were 
taken from the KPNO 4m FLAMINGOS catalog. In GOODS-S 
optical and near infrared observations were taken from the 
GOODS MUSIC catalog (Grazian et al. 2006; Santini et al. 
2009, UBVIzJHK). 



2.3. Removing AGNs 

To identify and remove X-ray AGN we use deep Chandra 
X-Ray observations, i.e., the 1 Ms maps for GOODS- 
S and the 2 Ms maps for GOODS-N (Alexander et al. 
2003). AGNs are identified as galaxies with either 
L x [0.5 - 8.0 keV] > 3 x 10 42 ergs -1 or a hardness ra- 
tio greater than 0.8 (ratio of the counts in the 2-8 keV to 0.5 - 2 
keV passbands) (Bauer et al. 2004). Even if it is well-known 
that AGN do also harbor star formation, we do not subtract the 
AGN contribution to the infrared light of those galaxies since 
such subtraction would be highly speculative at the present level 
of our knowledge. Instead we conservatively decide to remove 
all those galaxies from our sample. 



1 http://cdsweb.u-strasbg.fr/cgi-bin/qcat ?J/A+A/<volume>/<page> 



Fig. 2. Rest-frame 8 pm LF estimated at z ~ 2 using the 1/ V max method. 
Light blue circles and red squares represent the rest-frame 8 fim LF 
derived from our 24 fim sample with and without X-ray AGNs respec- 
tively. Empty triangles represent the rest-frame 8 fim LF derived by 
Caputi et al. (2007) at z ~ 2. Asterisks show the local reference taken 
from Huang et al. (2007) and the dotted line represents the best-fit to 
these data points with a double power law function with fixed slopes 
(i.e., (f> oc L~ ' 8 for L < L knee and <f> cc Zr 3 1 for L > L k „ ee ). The dashed 
line represents the best fit of the rest-frame 8 fim LF at z ~ 2 assum- 
ing that the shape of the rest-frame 8 fim LF remains the same since 
z ~ 0. The dark shaded area span all the solutions obtained with the x 1 
minimization method and compatible, within 1 <x, with our data points. 



2.4. The final infrared galaxy sample 

To construct our final infrared sample we first cross-match 
the 24 pm catalog with the X-ray observations. We find that 
~ 6%(13%) of the 24(70) pm sources contain an X-ray AGN. 
All those sources are excluded from our final infrared sample. 
Remaining sources are then matched with our spectroscopic and 
photometric redshift catalogs, using a matching radius of 1.5 " 
(i.e. ~ FWHM of the IRAC 3.6 pm observations). In case of 
multiple associations we select the closest optical counterparts. 
In GOODS-N (GOODS-S) 94% (92%) of the 24 /mi sources 
brighter than 30 yuJy have a spectroscopic and/or a photometric 
redshift and 80% (76%) of these sources have been detected in 
the near infrared. 46% and 37% of our 24 pm sources have a 
spectroscopic redshift in GOODS-S and -N respectively. 

All the different steps described previously are listed in 
the Table 1 . These steps yield to a final infrared galaxy sample 
containing 2644 and 138 sources detected at 24 pm and 70 pm 
respectively. In our redshift range of interest (i.e. 1.3 < z < 2.3) 
our infrared galaxy sample contains 706 sources detected at 
24 pm and only 8 sources detected at 70 pm. In this redshift 
range, the fraction of sources with a spectroscopic redshift is 
~ 20%. 

In Figure 1 (left), we compare spectroscopic and pho- 
tometric redshifts of 1670 sources detected at 24 pm and with 
both kinds of redshifts (i.e. spectroscopic and photometric 
redshifts). In the 1.3 < z < 2.3 redshift range, accuracy of 
the photometric redshifts is o"a z /(i +Z ) = 0.14 and Az/(1 + z) 
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Table 1. Redshift catalog properties 



Field 


Area 


Nb sources 


X-ray AGN 


No X-ray AGN 


# spec-z a 


# phot-z a - b 


# spec-z and/or phot-z a 




(arcmin 2 ) 


24//m/70 J um 


24/im/70/im 


24//m/70/im 


24/^111/70^111 


24 J um/70//m 


24 J um/70/im 


GOODS-S 


89 


870/44 


64/7 


806/37 


371/32 


378/5 


749/37 








7%/16% 


93%/84% 


46%/86% 


46%/ 14% 


92%/ 100% 


GOODS-N 


194 


2151/119 


134/12 


2017/107 


747/72 


1148/29 


1895/101 








6%/10% 


94%/90% 


37%/67% 


57%/27% 


94%/94% 



a percentages noted in these columns refer to the number of no X-ray AGN 

b number of sources which have a photometric redshift but no spectroscopic redshift 



has a median value of -0.002. These redshift uncertainties 
result in infrared luminosity uncertainties when converting 
24 /mi flux densities into L/« using the Chary & Elbaz (2001) 
library (see Figure 1, right). Since the redshift uncertainties 
are not Gaussian, the infrared luminosity uncertainties are also 
not Gaussian: wings of the real distribution extend further 
away than in a Gaussian distribution. As a result, to study the 
real impact of redshift uncertainties on the inferred infrared 
LF, one needs to introduce the real redshift distribution into 
Monte-Carlo simulations instead of using a redshift distribution 
with a Gaussian statistic. Such Monte-Carlo simulations have 
been done and are discussed in section 4. 1 . 

To illustrate the impact on the inferred LF of the sub- 
traction of the X-ray AGNs, we compute at z ~ 2 the rest-frame 
8 /mi LF with and without X-ray AGNs. The choice of this 
particular redshift and wavelength is motivated by the fact that 
at z ~ 2 the 24 /mi observations correspond to the rest-frame 
7.8 /mi. Therefore the extrapolation that needs to be applied to 
compute the 8 /im luminosity is negligible, nearly independent 
of the SED library used, as well as independent of the nature of 
the source (i.e., AGN or star-forming galaxy). In Figure 2 we 
present the rest-frame 8 /mi LF derived from the 1 /V max analysis 
(see section 4.1 for a precise description of this method) and 
using the Chary & Elbaz (2001) library. 

We note that the LF derived with and without the AGN 
are in total agreement except for the brightest luminosity bin. 
Even for this last luminosity bin the difference between these 
two LFs is relatively small and is of the order of ~ 0.2 dex. As 
a result, we can conclude that our particular choice for handling 
AGN will not affect our final results much and will not be able 
to explain the large discrepancies that will arise in later sections 
when we compare our LFs with others that have appeared in the 
literature. 

Our rest-frame 8 //m LF is in excellent agreement with the 
one inferred by Caputi et al. (2007), confirming the consistency 
of our 24 /mi sample. We note that the LF derived in our study 
extends to fainter luminosities since we are using a 24 /mi 
catalog that is ~ 3 times deeper. 



3. The Total Infrared Bolometric Correction 

In the local Universe tight correlations have been found be- 
tween monochromatic (eg, Li5^ m , i24/im and ... etc) and 
total infrared luminosities (Li R = L [8-1000 /mi]) of galaxies 
(Chary & Elbaz 2001). Based on these correlations SED libraries 
have been developed and extensively used to estimate the to- 
tal infrared luminosity of galaxies (e.g., Chary & Elbaz 2001; 



Lagache et al. 2003; Dale & Helou 2002, hereafter CE01, LDP, 
and DH 2 respectively). However, it is not clear that these lo- 
cal templates are suitable to describe the spectral properties of 
distant galaxies (Papovich et al. 2007; Daddi et al. 2007b; Le 
Borgne et al. 2009). To study this issue, we decided to character- 
ize the observed 24 /mi vs 70 /mi correlation and compare it to 
the predictions of standard SED libraries. Since at high redshift 
(i.e. z > 1.3) most of the 24 /mi sources are undetected at 70 /mi, 
the characterization of the 24 //m vs 70 /mi correlation could 
only rely on mean 70 //m properties obtained through a stack- 
ing analysis. This method, which has been extensively used in 
the last few years, gives reliable estimates of the typical 70 /mi 
flux density of a given galaxy population even below the detec- 
tion limit of current 70 /mi observations (e.g. Dole et al. 2006; 
Papovich et al. 2007; Magnelli et al. 2009). 

We first divided our 24 /mi sample into two redshift bins, 
of 1.3 < z < 1.8 and 1.8 < z < 2.3. Then, inside these redshift 
bins, we separated these 24 /im sources per luminosity bins of 
0.5 dex. For each 24 //m luminosity bin we stacked on the resid- 
ual 70 /im image all sources with no 70 /mi detection. The pho- 
tometry of the stacked image was then measured using an aper- 
ture radius of 16", a background from within annuli of 40"and 
60"and an aperture correction factor of 1 .705 (as discussed in the 
Spitzer observer's manual). Finally, the mean 70 /im flux density 
(F™^"") for a given 24 /mi luminosity bin was computed follow- 
ing Equation 1 : 



bin 



lOfim 



n + m 



(i) 



where F™^* is the stacked 70 /mi flux density of all 24 /im 
sources within this luminosity bin and undetected at 70 /im 
(sample which contains m sources); F. m is the 70 /im flux 
density of the ;th 24 /im sources within this luminosity bin and 
detected at 70 /im (sample which contains n sources). This 
procedure was performed using sliding 24 /mi luminosity bins 
with steps of 0. 1 dex. While such small sliding steps introduce 
correlations between our staking results, it avoids problems 
one might introduce by arbitrarily choosing some particular 
luminosity bins. Moreover, we note that since our staking anal- 
ysis probes a dynamic range of ~ 1.5 dex along the L24/im/(i+z) 
axis, there are always three independent measurements to 
characterize the typical L24 lim /(i +Z )-L 1 o lim /(\+z) correlation. 

Results of this stacking analysis are shown with filled 



2 We note that the Dale & Helou library is originally parametrized 
using the IRAS far-infrared colors (i.e., R60/ 100). Nevertheless the one 
used here has been parametrized a posteriori with L IR using the local 
^60/100 vs L IR correlation (Soifer & Neugebauer 1991). 
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Fig. 3. The 70 vs. 24 fim correlations as revealed by the observations and our stacking analysis in the two redshift bins considered in this study. 
The empty and filled squares represent the 70 vs. 24 fim correlations observed for sources individually detected at 70 fim with photometric or 
spectroscopic redshift respectively. The red diamonds show the results obtained using our stacking analysis (see text). For clarity the error bars 
of our stacking analysis are shown only for two points. These error bars are computed using a standard bootstrap analysis and an estimate of 
background fluctuation using a stacking analysis at random positions (for more detail see Magnelli et al. 2009). The empty triangle shows the 
median correlation found in the sample of Murphy et al. (2009). The thin solid lines, the dashed lines and the triple-dots-dash lines represent the 
expected correlations for the CE01, the LDP and the DH libraries, respectively, at the lowest and the highest redshift of each redshift bin. The 
red solid line represents the inferred 24/70 /jm correlation derived using a smooth linear interpolation between red diamonds and extended at high 
luminosities using 24 fim sources individually detected at 70 fim. At the bottom of each plot we present the fraction of 24 fim sources that are 
individually detected at 70 fim as a function of the 24/(1 +z) fim luminosity. 



red diamonds in Figure 3. In both redshift bins, we find 
that none of the usual SED libraries can reproduce the ob- 
served correlation between L24//m/(i+z) and L70/an/(i+z> At high 
24yum/(l+z) luminosities, standard SED libraries predict a 
higher Ljoiim/(i+z) to LaAiaaH\+z) ratio than is actually observed. 
These low Z/70/jm/(i+z) to L24/jm/(i+z) ratio could be reproduced by 
standard SED libraries but would correspond to SED templates 
with very low intrinsic 24/mi/(l + z) luminosities. In other 
words, the actual normalization of standard SED libraries which 
predict the increase of the L 7 o^ m /(i +Z ) to Z.24/im/(i+z) ratio with 
increasing 24/im/(l+z) luminosity is wrong. Instead, we 
observed that at high redshift the Lvo^m/fi+z) to L24/an/(i+z) ratio 
does not strongly depend on the 24/mi/(l + z) luminosity and 
that sources with high 24 fim/ (I + z) luminosity (or equivalently 
high infrared luminosity) have a Ljo^m/(i+z) to i24^m/(i+z) ratio 
typical of sources in the local universe with low 24/im/(l + z) 
luminosity. 

We note that among the few 24 fim sources detected at 
70 fim (open and filled squares in figure 3) those with high 
24yum/(l + z) luminosities (i.e. L24/jm/(i+z) > 8 x 10 11 L ) con- 
firm the discrepancy inferred using the stacking analysis. On the 
contrary, those sources with low 24/im/(l + z) luminosities are 
significantly closer to the predictions of standard SED libraries. 
This behavior is of course driven by selection effects and those 
sources only represent the high-end tail of the dispersion of the 
^24/«n/(i+z) vs Lj0fim/(i+z) correlation inferred through stacking. 

We note that the discrepancies that we find between the 
observed LiAnm/(i+z) vs Z/70//m/(i+z) correlation and predictions 
from standard SED libraries is quite different than what we 
found previously for galaxies at z < 1.3 using a similar stacking 
analysis (Magnelli et al. 2009). There, the CE01 models 
provided a reasonably good fit to the observed LiAnm/(i+z) vs 



Ljonm/a+z) correlation. 

Our findings indicate that the total infrared luminosity 
of a high-redshift galaxy cannot be inferred simply using its 
24 yum flux density and any of the three standard SED libraries 
considered here. We note that this analysis is fully consistent 
with the results presented by Daddi et al. (2007b, see their figure 
8) and first results obtained using Herschel data (Nordon et al. 
2010). At z ~ 2 and at L 24/J m/(i + z) = 8 X 10 1() and 3.5 x 10 11 L Q , 
Daddi et al. (2007b) find that Lj0f, m /(i+z) would be overestimated 
by a factor ~ 2 and ~ 10 using the CE01 library respectively, 
while we find a factor 2.3 and 12 respectively. 

In Papovich et al. (2007) and in Murphy et al. (2009) similar 
discrepancies were found and explained in part as the result of 
an increase in the PAH emission at any given infrared luminosity 
in high-redshift galaxies. In both studies, the infrared luminosity 
of galaxies were then derived using colder SED templates (i.e., 
like local galaxies with lower infrared luminosities) which fit 
the observed correlation. 

An alternative explanation for the discrepancies between 
SED libraries and the observed correlations can be the presence 
of an obscured AGN. Indeed, as suggested by Daddi et al. 
(2007a), the observed 24 fim flux density of a galaxy located 
at z ~ 2 might be dominated by the hot dust continuum from 
an obscured AGN and hence it might not be a robust SFR 
indicator as inferred from its disagreement with radio stacking 
and extinction corrected UV estimates. In the same study, Daddi 
et al. (2007b) show that in this redshift range (1.5 < z < 2.5) the 
SFRs derived from radio, extinction corrected UV and 70 fim 
flux density agree well. This suggests that unlike 24 fim, the 
70 fim passband is a good SFR indicator even in high-redshift 
galaxies. Thus, we will use the observed Z^/jm/ti+z) - ^70^m/(i+z) 
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Fig. 4. Different infrared bolometric corrections applied to a z ~ 1.8 
galaxy. Black square represents the observed 24 /im flux density while 
the black star represents the 70 /im flux density predicted using the 
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correlation. Dotted line and the dashed-dot line 



represent the unsealed CE01 templates corresponding respectively to 
the observed 24 /an and the predicted 70 /im flux densities. The dashed 
line represents the scaled CE01 template which best fit the 24 and 70 /mi 
flux densities of this galaxy. Ljj|, L^j and ll" R give the infrared luminos- 
ity derived by integrating the dotted line, the dashed-dot line and the 
dashed line respectively. 



correlation to infer Z/zo^m/a+z) for each 24 fim source and 
use this in turn to derive its total infrared luminosity and 
star-formation rate. 

Based on these two explanations we decided to derive 
the total infrared luminosity of each galaxy using two different 
methods, (i) For each 24 /im source we deduce its 70 //m flux 
density using the L24imi/(i+z) - £70imi/(i+z) correlation. Then we 
choose the CE01 template whose redshifted color best matches 
the derived 24 to 70 /mi flux ratio and renormalized it to match 
the 24 //m flux density of the source (i.e., ignoring the intrinsic 
luminosity normalization of the CE01 library; Dashed line of 
Fig. 4). The total infrared luminosity of this galaxy (hereafter 
Lj£) is estimated by integrating the SED curve of this template. 
In the following, we will assume that the uncertainties on the 
infrared luminosity estimated using this method is of order 0.2 
dex as measured by Murphy et al. (2009). We note that the 
i70iim/(i+z)/^24i/m/(i+z) ratio is nearly constant over the whole 
i24/jm/(i+z) luminosity range and corresponds to that expected 
for a CE01 template with an intrinsic Log(Lm [L G ] ) ~ 10.8. 
Such SED templates exhibit strong PAH features, (ii) For 
each 24 //m source, a 70 /im flux density is derived using the 
Luum/a+z) - L 1Qlim/{1+z) correlation. L IR (hereafter L™) is then 
simply estimated using the CE01 template (i.e., keeping the 
intrinsic luminosity normalization) matching the derived value 

Of £/70imi/(l+z)- 

Figure 4 illustrates these two different bolometric correc- 
tions for a galaxy situated at z = 1.8. We note that fitting the 24 
and 70 /im measurements together, allowing the renormalization 
of the CE01 SED templates, or using the luminosity-normalized 
CE01 library to fit the 70 /im alone give nearly the same results 
(for all our 24 //m sample we find < Log(Lf^/L™) >~ 0.04 and 

cr[Log(L^ I L 7 ^)] ~ 0.05). Indeed, both techniques involve the 
use of SED templates with lower intrinsic infrared luminosity 



than from the 24 //m alone. We also note that using the 24 //m 
flux density alone and the luminosity-normalized CE01 library 
we would have overestimated the total infrared luminosity of 
this galaxy by ~ 0.4 dex. 



4. The Infrared Luminosity Function 

4.1. Methodology 

The infrared LFs are derived using the standard l/V max method 
(Schmidt 1968). The comoving volume associated with any 
source of a given luminosity is defined as V max = V zmax - V zmi „ 
where zmin is the lower limit of the redshift bin, and zmax is 
the maximum redshift at which the object could be seen given 
the flux density limit of the sample, with a maximum value 
corresponding to the upper limit of the redshift bin. For each 
luminosity bin, the LF is then given by 



-7 

AL^ 



1 



V max ,i X Wi 



(2) 



where V max is the comoving volume over which the /th galaxy 
could be observed, AL is the size of the luminosity bin, and w, 
is the completeness correction factor of the /th galaxy, w, equals 
1 for sources brighter than S241™ ~ 100/Jy and decreases at 
fainter flux densities due to the incompleteness of the 24 /im 
catalog. These completeness correction factors are robustly 
determined using Monte-carlo simulations (Magnelli et al. 
2009) and reach at most a value of 0.8. None of the conclusions 
presented here strongly rely on this correction. 

Uncertainties in the infrared LF values depend on 
photometric redshift uncertainties (see Figure 1). In particular, 
catastrophic redshift errors (i.e. Az/(1 + z) > 0.15), which can 
shift a low redshift galaxy to higher redshift and vice versa, can 
modify the number density of LIRGs and ULIRGs in a given 
redshift bin. To estimate the effect of these catastrophic redshift 
errors on the derived infrared LF, one needs to have a complete 
census on the population of infrared galaxies at all redshifts. To 
simulate this, we first generate a reference catalog (i.e., an ideal 
sample with no redshift uncertainties), then from this reference 
catalog we generate 1000 mock catalogs with realistic redshift 
uncertainties and finally we compare the infrared LF inferred 
from the reference catalog with the infrared LF inferred from 
the 1000 mock catalogs. The key point of these simulations is 
to introduce redshift uncertainties which accurately reproduce 
the observed distribution of spectroscopic versus photometric 
redshifts (see Figure 1) instead of using a standard Gaussian 
distribution which would not be realistic (see discussion in 
section 2.2). 

The reference catalog is constructed as follows. We 
first start from a simulated catalog generated from the model of 
Le Borgne et al. (2009) that best fits number counts of sources 
at 15, 24, 70, 160, 850 /an. This catalog, which contains all 
infrared sources that should be observed over a field of 283 
arcmin 2 with < z < 5, reproduces the observed infrared LF 
up to z ~ 1.3 (e.g., see figure 12 of Magnelli et al. 2009). We 
only keep sources with 0<z<1.3orz> 2.3. Then from 
our observed infrared LF (section 4.2) we construct and add 
to this catalog all infrared sources that should be observed at 
1.3 < z < 2.3 over a field of 283 arcmin 2 . This catalog, which 
by construction reproduces all the observed infrared LF from 
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Fig. 5. Results of our Monte Carlo simulations. Black lines represent the infrared LF that one would have inferred using an ideal sample with 
no redshift or bolometric correction uncertainties. Red squares show the mean infrared LF inferred using our 1000 mock catalogs. Error bars 
correspond to the dispersion observed in our 1000 Monte Carlo simulations. 
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Fig. 6. Infrared LF estimated in two redshift bins with the 1/V„ M _ V method. Red and dark blue squares are obtained using L t ' R and LjjJ respectively 
(see text). Asterisks show the local reference taken from Sanders et al. (2003) and the dotted line represents the best-fit to these data points 
with a double power law function with fixed slopes (see text). The dashed line represents the best fit of the infrared LF assuming that the shape 
of the infrared LF remains the same since z ~ 0. The dark shaded area span all the solutions obtained with the x 2 minimization method and 
compatible, within 1 a, with our data points. In the first redshift panel, we reproduce in green, blue, yellow and red the best fit of the LF obtained 
at 0.4 < z < 0.7, 0.7 < z < 1.0, 1.0 < z < 1.3 (Magnelli et al. 2009), and 1.8 < z < 2.3 respectively. 



z = up to z ~ 2.3, will be our reference catalog. 

Starting from this reference catalog we create 1000 mock 
catalogs which contain the same number of sources as the 
original one but we attribute to each source a new redshift 
randomly selected to reproduce the observed distribution of 
spectroscopic versus photometric redshifts (see Figure 1). To 
take into account the bolometric correction uncertainties not 
associated with the photometric redshift errors, we attribute 
to each source a new infrared luminosity selected inside a 



Gaussian distribution centered at the original source luminosity 
and with a dispersion of 0.2 dex (see Section 3). Using these 
1000 mock catalogs, we then compute the infrared LF and study 
the difference between these infrared LFs and the infrared LF 
derived from the reference catalog. 

Using these Monte Carlo simulations, we find only 
small systematic offsets between the real infrared LF and the 
one inferred in presence of redshift uncertainties (see figure 
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Fig. 7. Evolution of (pknee and L^„ Cf , as function of redshift. Data points 
below z = 1.3 are taken from Magnelli et al. (2009). 



5). At z ~ 2 and at faint luminosities, we underestimate the 
LF values by at most 0.1-0.15 dex while at bright luminosities 
we overestimate the LF values by at most 0.1-0.2 dex. These 
systematic offsets are smaller than the total uncertainty in each 
luminosity bin (~ 0.25 dex) defined as the quadratic sum of 
the Poissonian error (oc 1/ -yW) and the dispersion given by the 
Monte Carlo simulations. As a result, in the following we do not 
correct the inferred LF for these systematic offsets. 

4.2. Results 

In figure 6, we present the infrared LF derived in two redshift 
bins (1.3 < z < 1.8 and 1.8 < z < 2.3) using our two different 
infrared bolometric corrections (red and blue squares for 
and L 1 ^ respectively, see Table B.6 and B.7). First, we note that 
the infrared LFs derived using these two different bolometric 
corrections are in very good agreement and certainly within the 
error bars. Here on, we will refer to the LF derived using Lj' 
as the infrared LF. Indeed we have a better understanding of the 
uncertainties of this technique and this bolometric correction 
was found to be reliable up to z ~ 2 by Murphy et al. (2009). 

We take as a local reference the infrared LF derived by 
Sanders et al. (2003), showing their data points (stars) and their 
best double power law fit (i.e., <p oc Zr° 6 for Log(L/L G )< 10.5 
and <p oc Zr 2 - 15 for Log(L/L )> 10.5). Then, using a y 1 
minimization, we fit our infrared LFs with the same function, 
fixing the power law slopes at their z ~ values and leaving 
Lk„ ee and fyknee as free parameters (see Table 2). The shaded 
regions present all the solutions which are compatible with the 
data within 1 cr. These shaded regions extend to luminosities 
lower than our current observations and strongly depend on the 
assumption we made on the shape of the infrared LF, i.e. that 
it remains the same since z ~ 0. We will see latter on that this 
assumption is quite consistent with observational constraints 
obtained on the low luminosity end of the infrared LF (Reddy 
et al. 2008). 

The evolution of Lk nee and <pk nee between z = 1.3 and 
z = 2.3 is compared in Figure 7 with the evolution found at 



Fig. 8. The infrared LF at z ~ 2 obtained in this work (dark shaded 
area and dashed line) as compared with the determinations of other 
authors. The infrared LF at z ~ 2 obtained by Caputi et al. (2007, we are 
using their double exponential function) is represented by the dashed 
dotted line. Empty circles represent the infrared LF at z ~ 2 inferred by 
Perez-Gonzalez et al. (2005). Empty diamonds represent the infrared 
LF at z ~ 2.3 inferred by Reddy et al. (2008). Empty triangles represent 
the infrared LF at z ~ 2.5 inferred by Chapman et al. (2005). Filled 
black diamonds and empty squares represent the infrared LF that we 
would have inferred using our 24 /jm sample and the unsealed LDP or 
CE01 SED libraries respectively. The horizontal dashed line presents 
the source density below which the number of sources in the volume of 
GOODS and in a luminosity bin of 0.5 dex is less than 2. The vertical 
dashed line represents the corresponding luminosity using the best fit of 
our infrared LF. 

lower redshifts by Magnelli et al. (2009). Assuming that the 
shape of the LF remains the same since z ~ 0, we express 
the evolution of the infrared LF as p(L,z) = g(z)p(L/ f(z),Q), 
where g(z) and f(z) describe the density and the luminosity 
evolution through g(z) — (1 + z) p and f(z) = (1 + z) q . Between 
z — and z ~ 1 the redshift evolution consists mainly of a 
slight density evolution proportional to (1 + z )-°- 9±0 - 6 and a 
luminosity evolution proportional to (1 +z) 3 5±0 5 . Then between 
Z ~ 1 and z ~ 2 we observe a density evolution proportional to 
(1 +z) _1 - 1±L5 associated with a luminosity evolution proportional 
to (1 + z ) 10±0 - 9 . In comparison, at z > 1 Caputi et al. (2007) 
find a density evolution proportional to (1 + z )- 3 - 9±l ° and a 
luminosity evolution proportional to (1 + f) 2 - 2±0 - 5 . The evolution 
of the infrared LF that we find at z > 1.3 is more gradual than 
that derived by Caputi et al. (2007) and is nearly consistent with 
no evolution. 

Appendix A presents the evolution of the rest-frame 
8, 15, 25, 35 //m LFs. Rest-frame luminosities of each source 
are derived using the SED that we used to compute its L^. 



4.3. Comparison with previous work 

In Figure 8 we compare our results at z ~ 2 with the infrared LF 
inferred in various previous studies. There is a clear disagree- 
ment between our results and the infrared LF derived by Caputi 
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Rprkhiff 

JXCUMlllL 




"2 












Log(L ) 




z~0 


-0.60 


-2.20 


10.48 ± 0.02 


-2.52 ± 0.03 


0.4<z<0.7 b 


-0.60 


-2.20 


11.19 ±0.04 


-2.84 ± 0.06 


0.7 < z < 1.0 b 


-0.60 


-2.20 


11.37 ±0.03 


-2.65 ± 0.05 


1.0<z< 1.3 b 


-0.60 


-2.20 


11.69 ±0.06 


-2.91 ±0.10 


1.3 <z< 1.8 


-0.60 


-2.20 


11.84 ± 0.13 


-3.07 ± 0.20 


1.8<z<2.3 


-0.60 


-2.20 


11.83 ± 0.13 


-3.04 ± 0.22 



a Fixed slopes 

b These parameter values are taken from Magnelli et al. (2009) 



et al. (2007). This discrepancy arises from the fact that the 
bolometric corrections used in Caputi et al. (2007, i.e. with the 
LDP library) do not take into account the SED evolution that we 
observe at high redshift. Using the same bolometric corrections 
as Caputi et al. (2007) does lead us to similar results (see black 
diamonds of Figure 8). The disagreement of our results with the 
LF from Perez-Gonzalez et al. (2005) is even larger because 
they compute their bolometric corrections using the CE01 
library (see open squares in Figure 8), i.e. the standard SED 
library which exhibit the largest discrepancies with the observed 
^/jm/d+z) - £7(W(i+z) correlation at z > 1.3 (see Section 3). 
Rodighiero et al. (2010) have also derived the z ~ 2 infrared 
LF using deep 24 pm observations of the VVDS-SWIRE 
(S24/,m > 400 //Jy) and GOODS (S^m > 80 pJy) fields. Data 
points from this study are not shown in Figure 8 since they are 
very similar to that from Caputi et al. (2007) (see figure 15 of 
Rodighiero et al. 2010). 

In Figure 8, we also compare our results with the infrared 
LF inferred at z ~ 2.3 by Reddy et al. (2008) using observations 
of UV-selected star-forming galaxies. This study derived SFR 
and dust reddening from the UV rest-frame observations 
calibrated by comparison to 24 pm photometry for brighter 
sources. The UV-derived extinction was used to compute the 
expected infrared emission from galaxies fainter than the 24 pm 
detection limit and hence provide an extension of the IR LF to 
fainter luminosities. We find good agreement between our LF 
and that of Reddy et al. in the range of luminosities where the 
two studies overlap. At faint luminosities (log(L/ R /L ) < 11) 
our best fit LF falls somewhat below that derived by Reddy 
et al. (2008), although they are still consistent within the 
uncertainties. We have no direct measurements with Spitzer at 
such faint luminosities and rely upon an extrapolation based on 
a faint-end slope fixed at its z ~ value. As detail in section 4.4 
this disagreement at faint luminosities has nearly no impact on 
the SFR density inferred at z ~ 2. Indeed, the integrated SFR 
density of the universe at z ~ 2 computed from our LF agrees 
with that derived by Reddy et al. (2008). 

Finally we compare our results with the infrared LF inferred 
at z ~ 2.5 by Chapman et al. (2005) using submillimeter obser- 
vations 3 . The luminosity range probed by Chapman et al. (2005) 
is not constrained by our study since the comoving volume 
probed by the GOODS survey is too small (i.e., fewer than 2 



3 S. Chapman (private communication) confirms that the far-infrared 
luminosities for the submm galaxy LF reported in Table 6 of Chapman 
et al. (2005) are integrated over the wavelength range 8- 1 100 jim, nearly 
the same as the range 8-1000 jim that we adopt here. He also verifies 
that the luminosity column of that table should include an unspecified 
factor of h 2 . We have converted the data from Chapman et al. (2005) 
Table 6 to the value H = 70 km s -1 Mpc~ 3 that we adopt in this paper. 



sources would be present in this volume for a luminosity bin of 
ALog(L IR ) = 0.5 dex). We note however that the extrapolation 
of our infrared LF to high luminosities is consistent with the 
estimates of Chapman et al. (2005). 



4.4. Discussion 

We derive the evolution of the comoving number density 
of LIRGs and ULIRGs by integrating the infrared LF at 
z = 1.55 + 0.25 and z = 2.05 ± 0.25. We then combine these 
estimates with the evolution found at < z < 1 .3 by Magnelli 
et al. (2009) (Figure 9 left). We find that the number densities 
of LIRGs and ULIRGs between z ~ 1.3 and z ~ 2.3 are 
nearly constant. The number density of ULIRGs at z ~ 2 
(7-5^40 x 10" 5 Mpc~ 3 ) agrees with estimates made by Daddi 
et al. (2007b) (~ 10 x 10~ 5 Mpc~ 3 ) using UV observations 
calibrated against radio and other non-24 pm data. However, if 
we compare our number density with the Daddi et al. (2007b) 
estimate obtained combining UV and 24 pm observations 
(16^g 2 x 10~ 5 Mpc~ 3 ) we find a clear disagreement. 

Figure 9 (right) presents the evolution of the comoving 
infrared luminosity density (IR LD; or, equivalently, SFR 
density under the assumption that the SFR and L IR are related 
by Eq. 3 for a Salpeter IMF, i.e. 0(m)oc m~ 2 35 , between 0.1-100 
M Q , Kennicutt 1998) produced by ULIRGs, LIRGs, and by 
galaxies with Lm < 10 n L o (hereafter called "normal" galaxies 
by analogy with ordinary spiral galaxies at z ~ 0, although we 
note that at high redshift, LIRGs and ULIRGs themselves are 
sufficiently common to be considered "normal" for their epoch). 
We find a slight decrease of the IR LD from z = 1.3 to z = 2.3 
due to a decrease in the contribution of LIRGs and "normal" 
galaxies. At z ~ 2 the IR LD of the Universe is still dominated 
by LIRGs and not by ULIRGs, contrary to previous claims (e.g., 
Perez-Gonzalez et al. 2005). Using the best fit of our infrared 
LF, we infer that at z ~ 2 LIRGs and ULIRGs have an IR LD of 
4.5 x 1O 8 L Mpc~' and 1.5 x 1O 8 L Mpc~' respectively and that 
they account for 49% and 17% of the total IR LD respectively. 

SFR [M Q yr 1 ] = 1.72 x 10- 10 L IR [L ] (3) 

We compare our estimates with the evolution derived by Caputi 
et al. (2007) (Figure 10 left). At z ~ 2 we find that the IR LD of 
ULIRGs estimated from our best fit is a factor of ~ 1.8 lower 
than that estimated by Caputi et al. (2007). For LIRGs, for 
which we still have a good constraint due to our deep 24 pm 
sample, our estimate of their IR LD is a factor of ~ 1.5 higher 
than that inferred by Caputi et al. (2007). Nevertheless, we note 
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Redshift Redshift 

Fig.9. (Left) Evolution up to z ~ 2.3 of the comoving number density of "normal" galaxies (i.e. 10 7 L G < Lj r < 10" L G ; black filled triangles), 
LIRGs (orange filled diamonds) and ULIRGs (red filled stars). The green circles represent the total number of galaxies which are above the 24 pm 
detection limit of the surveys presented here, i.e. L ir > L?"* llmlt . The z ~ points are taken from Sanders et al. (2003). (Right) Evolution of the 
comoving IR energy density up to z ~ 2.3 (upper striped area) and the relative contribution of "normal" galaxies (yellow filled area), LIRGs 
(orange filled area) and ULIRGs (red filled area). The areas are denned using all the solutions compatible within 1 a with the infrared LF. Black 
arrows show the comoving IR energy density derived by stacking the 70 pm image at all IRAC sources positions (S'^^ m > 0.6 /Jy). The axis on 
the right side of the diagram shows the evolution of the SFR density under the assumption that the SFR and Li R are related by Eq. 3 for a Salpeter 
IMF. 
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Fig. 10. (Left) Evolution of the comoving IR energy density up to z ~ 2.3. Blue empty circles represent the results obtained by Caputi et al. (2007) 
for the global evolution of the comoving energy density (solid line) and the relative contribution of "normal" galaxies (dot line), LIRGs (dashed 
line) and ULIRGs (dot dashed line). Filled black star represents the comoving IR energy density of the Universe inferred at z ~ 2.3 by Reddy et al. 
(2008) while open star shows the relative contribution of LIRGs. Filled areas are as in Figure 9. (Right) Evolution of the comoving SFR density 
up to z ~ 2.3 assuming that SFR and Lir are related by Eq. 3 for a Salpeter IMF. Filled areas are as in Figure 9. The dotted line represents the 
SFR measured using the UV light not corrected for dust extinction (Tresse et al. 2007). The dashed line represents the total SFR density defined 
as the sum of the SFR density estimated using our infrared observations and the SFR density obtained from the UV light uncorrected for dust 
extinction. Light blue diamonds are taken from Hopkins & Beacom (2006) and represent the SFR densities estimates using various estimators. 
Dark blue triangles represent the SFR density estimated by Seymour et al. (2008) using deep radio observations. Green circles represent the SFR 
density estimated by Smolcic et al. (2009) using deep 20 cm observations and dark blue squares represent the relative contribution of ULIRGs to 
this SFR density. 



that if we take into account the range of IR LD defined by all 
solutions compatible within l<x with our data points, our IR 
LDs of LIRGs and ULIRGs are compatible with estimates from 
Caputi et al. (2007). On the other hand, the IR LD estimated by 
Caputi et al. (2007) for galaxies with Lir < 10 11 L Q is far below 
our estimate since they used a flatter faint-end slope for their 
infrared LF. We believe that our estimate is more reliable since 
we are using a 24 /mi catalog that is ~ 3 times deeper than that 
used by Caputi et al. (2007). We also note that the extrapolation 
of our infrared LF to these faint luminosities is corroborated by 



the infrared LF inferred by Reddy et al. (2008). 

We also compare ourz ~ 2.05 IR LD values with the z ~ 2.3 
estimates of Reddy et al. (2008). We note that this comparison 
is not straightforward because their IR LD needs to be slightly 
corrected prior to be compared with our work. Indeed, since 
they cannot constrain with their sample the contribution of 
ULIRGs to the IR LD, they use the value derived by Caputi 
et al. (2007). By replacing the Caputi et al. (2007) estimates 
by our value we compute the correct IR LD of Reddy et al. 
(2008), i.e., 10.0 + 0.2 x 10 8 L Mpc" 3 . We find that the IR LDs 



Magnelli et al.: Infrared Luminosity Density at < z < 2.3 



11 



derived by Reddy et al. (2008) for LIRGs and for the integrated 
population of all galaxies agree with our values. This agreement 
reinforces the idea that at z ~ 2 UV corrected for extinction is a 
good SFR indicator (Daddi et al. 2007b). 

In order to get a complete census on the SFR history 
we need to take into account the contribution of unobscured 
UV light. The unobscured SFR density {dotted line in Figure 
10 right) is taken from Tresse et al. (2007) and corresponds to 
the SFR density inferred using UV observations not corrected 
for extinction. The total SFR density {dashed line in Figure 10 
right) is then defined as the sum of the unobscured SFR density 
traced by the direct UV light and the dusty SFR density traced 
by the infrared emission. We find that the relative contribution 
of unobscured UV light to the cosmic SFR density evolves 
nearly in parallel with the total one and accounts for ~20% 
of the total SFR density. Globally, the cosmic star-formation 
history that we derived is consistent with the combination of 
indicators, either unobscured or corrected for dust extinction, 
as compiled by Hopkins & Beacom (2006). We also notice a 
very good agreement between the cosmic star-formation history 
derived in our work and the ones derived by Seymour et al. 
(2008) using deep radio observations. 

Our estimates can suffer from several uncertainties. 
Especially the contribution of "normal" galaxies to the IR LD 
comes from the extrapolation to low luminosities of the infrared 
LF where we have no constraints. To cross check our results 
we compute a lower limit on the comoving IR LD by stacking 
70 pm images at the positions of all IRAC sources in each 
redshift bin of interest (i.e., S' 3 R ^ m > 0.6//Jy; up arrows in 
Figure 9 right). This analysis is possible because the correlation 
between L7o^m/(i+z) and L m is quasi-linear at this redshift, and 
hence 2S(70//m) oc I.L IR . The stacking result is fully consistent 
with the value based on the integration of the extrapolated best 
fit to our infrared LF. 

As discussed in Section 3, the role of obscured AGN 
on the estimate of the infrared LF is still uncertain. Such results 
will be debated until the Herschel infrared space observatory 
provides accurate far-infrared measurements for faint, high- 
redshift galaxies. However, as shown by Murphy et al. (2009) 
using IRS spectroscopy, the mid-IR continuum from an AGN 
appears to scale with increasing 24 pm luminosity. As a result, 
the removal of any additional contribution from obscured AGN 
activity will only steepen the bright-end of the infrared LF. 
This would reinforce our main result which is the fact that at 
z ~ 2 previous studies have overestimated the number density 
ofULIRGs. 



5. Conclusion 

For the first time we take advantage of deep far-infrared obser- 
vations to derive the evolution of the infrared luminosity density 
(or equivalently SFR density under the assumption that the SFR 
and Lir are related by Eq. 3 for a Salpeter IMF) over the last 
4/5ths of the cosmic time (see Figure 9 right). Using the deepest 
24 pm (~ 3 times deeper than any previous studies) and 70 /im 
observations made by Spitzer and a careful stacking analysis we 
are able to calibrate a new infrared bolometric correction based 
on the renormalization of local SED templates reproducing the 
observed L24/(i +Z ) vs L 7 o/(i +Z ) correlations. This new bolometric 



correction is a key result of our analysis, since previous studies 
on the infrared LF at z ~ 2 (e.g., Perez-Gonzalez et al. 2005; 
Caputi et al. 2007) did not account for the evolution of the in- 
frared SEDs that we observe in high-redshift galaxies. Our main 
result is that we find a flattening of the SFR density between 
z = 2.3 and 1.3, and that the comoving number density of LIRGs 
and ULIRGs remain nearly constant over this redshift range. At 
z ~ 2 the SFR density is still dominated by LIRGs and not by 
ULIRGs contrary to previous claims. The flattening of the SFR 
density at z > 1 reinforces the idea that at this redshift we ob- 
serve a change of the properties of star-forming galaxies as al- 
ready suggested by the reversal of the star formation-density re- 
lation at z ~ 1 (Elbaz et al. 2007). 

The evolution of the SFR density of the Universe provides 
a strong constraint on the main mechanism which triggers the 
SFR in galaxies. At z < 1 the decrease of the SFR density might 
be driven by a gradual gas exhaustion as suggested by the con- 
tinuous decrease of the SFR vs M* relation in this redshift range 
(Noeske et al. 2007; Elbaz et al. 2007). Between z ~ 2 and z ~ 1, 
the relatively constant SFR density still needs to be understood 
in the framework of large-scale structure formation, merging, 
and/or AGN activity. 

The main limitation of our study is the uncertainty on the in- 
fluence of obscured AGN on the infrared bolometric correction 
to be applied to bright 24 pm sources. This influence will soon 
be assessed using new far-infrared observations from Herschel. 
In particular, the GOODS-Herschel Open Time Key Programme 
(PI: David Elbaz) will reach the faintest flux limits at 100 pm in 
an ultradeep field within GOODS-S, expected to provide indi- 
vidual measurements for most z ~ 2 galaxies detected at 24 pm 
by Spitzer, where here we could only derive average values 
based on 70 pm stacking. This should help disentangle the con- 
tributions of AGN and star formation for sources over a broad 
swath of the high-redshift infrared luminosity function. 
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Appendix A: The Rest-frame 8 jum, 15 yum, 25 yum 
and 35 yum LFs 

We aim to derive the rest-frame 8 /im, 15 fim, 25 fim and 35 pm LF from our 
24 fim sample. This is done here for several reasons. While the derivation of 
the rest-frame 8 fim LF has already been addressed in some previous studies 
(Caputi et al. 2007), our 24 fim sample reaches flux limits ~ 3 times fainter, 
providing improved constraints on the LF break. The rest-frame 15 pm LF 
provides continuity with what we have computed in Magnelli et al. (2009). 
The rest-frame 25 /im LF was not derived in Magnelli et al. (2009) but it does 
have several points of interest. First, it reduces the i-correction that one has to 
apply when observing infrared sources at z ~ 2 using a 70 /im passband or at 
z ~ 3 using a 100 fim passband. Hence, this rest-frame LF could be compared 
with future 70 fim and 100 /im observations made using the Photodetector 
Array Camera and Spectrometer (PACS) instrument onboard the Herschel 
satellite. Second, the rest-frame 25 /im LF offers a means to compare directly 
with IRAS 25 /im observations of galaxies in the local universe. Third, and 
perhaps most importantly, the rest-frame 25 /im luminosity has been shown to 
provide a reliable and fairly direct measurement of star formation in galaxies 
(e.g., Calzetti et al. 2007). Hence, one may want to use these measurements 
in the future to directly derive the SFR distribution function without using 
intermediate bolometric corrections. The rest frame 35 /im LF also provides 
continuity with quantities derived in Magnelli et al. (2009) and corresponds 
at z, ~ 2 to the observed 100 /im flux density. As a result, this rest-frame LF 
would be a standard comparison for the z ~ 2 LF computed using PACS 100 /im 
observations. 

The various rest-frame luminosities are derived using the same method 
as that used to compute ll" R . For each 24 /im source we deduce its 70 /im 
flux density using the L24pm/(l+z) - £70//m/(l+z) correlation. Then we choose 
in the CE01 library the scaled template which best fits these 24 and 70 /im 
fluxes densities. Using this scaled template we then compute the rest-frame 
luminosities of this galaxy in our four passbands of interest (i.e., at 8 fim, 15 
/im, 25 fim and 35 /im). The corresponding rest-frame LFs are computed using 
the 1 IV m ax method. All the LF are then fitted using a double power law function 
with fixed slopes estimated using a bivariate method (see Table B.l). 

Figure B.l presents the rest-frame 8 /im LF derived in our two redshift 
bins using the l/V max method. These LF are compared with the local reference 
of Huang et al. (2007) and with the z ~ 2 LF derived by Caputi et al. (2007). 
Our estimates agree well with those of Caputi et al. (2007), but extend to lower 
luminosities. 

Figure B.2 presents the rest-frame 15 /im LF derived in our two redshift 
bins using the 1 /V 'max method. Comparing to the local reference of Xu (2000) 
and with the z, ~ 0.55, z ~ 0.85 and z ~ 1.15 reference of Magnelli etal. (2009), 
we note the strong evolution of this LF with redshift. 

Figure B.3 presents the rest-frame 25 /im LF derived in our two redshift 
bins using the l/V max method as well as the local reference of Shupe et al. 
(1998). 

Finally Figure B.4 presents the rest-frame 35 /im LF derived in our two 
redshift bins using the \IV max method as well as the local reference derived 
from Shupe et al. (1998), and the LFs derived in Magnelli et al. (2009) at 
z ~ 0.55, z ~ 0.85 and z ~ 1.15. 



Appendix B: Source list 

At the resolution of Spitzer most of the sources in our fields are point sources 
(i.e FWHM~ 5.9" and 18"at 24 fim and 70 /im respectively). Therefore, to 
derive their photometry we decided to use a PSF fitting technique that take into 
account, as prior information, the expected position of the sources. Starting 
from IRAC positions (GOODS-N: GOODS legacy program, Dickinson et al., 
in preparation; GOODS-S: SIMPLE catalog, Damen et al. 2011) we extract all 
24 /im sources. Then, using our 24 /im catalogs, we extract all 70 /im sources. 
This method deals with a large part of the blending issues encountered in dense 
fields and allows straightforward multi-wavelength association between near-, 
mid- and far-infrared sources. The disadvantage of this method is that we have 
to assume that all sources present in our mid-infrared images have already been 
detected at IRAC wavelengths. In our case this assumption is true because our 
IRAC 3.6 /im data are 30 times deeper than our current 24 fim observations and 
that the typical S 24^1/^3 6^ m ratio spans the range [2-20]. 

In this online material we release our complete 24 fim and 70 fim 
source catalogs for both GOODS fields. These catalogs expend below the 80% 
completeness limit, and cover the full area (approximately 10' X 16') of the 
GOODS-S region (i.e., not only the smaller 10' X 10' region with the deepest 
70 /im imaging that is used for the analysis in this paper). The noise level in 
the GOODS 24 /im data is homogeneous over most of the field, with some 
degradation near the edge where the exposure time is somewhat reduced. We 



restrict our release to regions with exposure time higher than 9500 s per pixel. 
This limit corresponds to a quarter of the exposure time of the deep inner region 
(~ 38000 s). This degradation does not really affect the depth of our 24 /im 
catalogs in that region since uncertainties are still dominated by confusion 
(cr™ p e « 2 x cr™2'„ ~ 6/Jy < o- co "/""°" » 7/dy). At 70 /im, the noise level 
is roughly uniform throughout GOODS-N (~ 12000 s per pixel). However, in 
GOODS-S the deepest 70 /im data, with noise similar to those in GOODS-N, 
are limited to a region approximately 10' X 10' in extent. The outer region 
portions of the GOODS-S field have somewhat shallower 70 /im data (~ 6 000 s 
per pixel). 

At 24 fim, sources are detected using an empirical 24 fim PSF con- 
structed with isolated point like objects present in the mosaic. At 70 fim no 
reliable empirical PSF could be constructed because only a few isolated sources 
could be found in each map. We then decided to use the appropriate 70 fim Point 
Response Function (PRF) estimated on the extragalactic First Look Survey 
mosaic (xFLS; Frayer et al. 2006a) and available on the Spitzer web site. At both 
wavelengths an aperture correction is applied to all flux densities to account for 
the finite size of our PSFs. Those aperture corrections are taken from the Spitzer 
data handbook. 

Calibration factors used to generate the final 24 and 70 fim mosaics are 
derived from stars, whose SED at these wavelengths are generally very different 
from those of distant galaxies. Hence, color-corrections have to be applied to 
all flux densities (at most ~ 10%). In the catalogs released here, 70 fim flux 
densities have been color-corrected using a systematic and standard correction 
of 1.09 (see Spitzer data handbook). This 70 fim color-correction is computed 
for distant galaxies with dust temperature of ~ 40 K. This color-correction 
differs from those applied in our study and which take into account the redshift 
of each source (see discussion in Section 2.1). No color-corrections are applied 
to our 24 fim flux densities since, for those data, color-corrections are more 
strongly dependent on the redshift of the source. Indeed, 24 fim data probes 
different part of galaxy SED as function of the redshift (black-body emission of 
dust or PAH emission). 

Our 24 fim and 70 fim data are the deepest observations taken by 
Spitzer and have been designed to reach the confusion limit of this satellite. 
Flux uncertainties are therefore a complex combination of photon and confusion 
noise. In order to estimate these complex flux uncertainties and to characterize 
the quality of our 24 /im and 70 fim catalogs we use two different approaches. 
First, we compute the noise of each detection using our residual maps. Second, 
we estimate a statistical flux uncertainty based on extensive Monte-Carlo 
simulations. 

Noises estimated on residual maps correspond to the pixel disper- 
sion, around a given source, of the residual map convolved with the appropriate 
PSF. This method has the advantage of taking into account the rms of the map 
and the quality of our fitting procedure. These noise estimates are given in our 
released catalogs as <j map . These estimates are almost equal to the rms of our 
maps, i.e., cr map ~ 3 /ily/beam at 24 fim in both GOODS-N and GOODS-S, 
0.3 mly/beam at 70 fim in GOODS-N and in the deepest region of GOODS-S, 
and 0.45 mly/beam at 70 fim in the shallowest region of GOODS-S. 

In order to estimate the effect of confusion noise we performed ex- 
tensive Monte-Carlo simulations. We added artificial sources in the 24 fim 
and 70 fim images with a flux distribution matching approximately the measured 
number counts (see Frayer et al. 2006b; Papovich et al. 2004). To preserve the 
original statistics of the image (especially the crowding properties) the numbers 
of artificial objects added in the image was kept small (we only added 40 sources 
into the 24 fim images and 4 sources into 70 fim images). We then performed 
our source extraction method and compared the resulting photometry to the 
input values. To increase the statistic, we used repeatedly the same procedure 
with different positions in the same field. For each field we introduced a total of 
20000 artificial objects. Results of these Monte-Carlo simulations are shown in 
Figure 1 of Magnelli et al. (2009) and are summarized hereafter. 

From these Monte-Carlo simulations we derive three important quantities: 
the photometric accuracy, the completeness and the contamination of our 
catalogs as function of flux density. Completeness is define as the fraction 
of simulated sources extracted with a flux accuracy better than 50%. The 
contamination is defined as the fraction of simulated sources introduced with 
S < 2cr" mp which are extracted with S > 3o Jnap . 

Using these Monte-Carlo simulations, we find that in both GOODS fields 
our 24 fim catalogs are 80% complete at ~ 30 /ily. At this flux density, the flux 
accuracy is better that 20% and the contamination is ~ 10%. The flux accuracy 
of our source extraction reaches 33% around 20 /ily. This limit could be defined 
as the "real" 3cr" m " limit of our data because this estimate take into account 
confusion noise. At 20 /ily, the completeness of our catalog is ~ 40% and the 
contamination is ~ 15%. 

For our deep 70 fim data in GOODS-N and -S, Monte-Carlo simulations 
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show that our catalogs are 80% complete at 2.5 mJy. The 33% flux accuracy 
is reached at 2 mJy with a completeness of ~ 50% and a contamination of 
~ 15%. For the shallow 70 fim data of GOODS-S, the 80% completeness limit 
is reached at 3 mJy, and the 33% flux accuracy is reached at 2.5 mJy. At 2.5 
mJy, the completeness is 45% and the contamination is 15%. 

Flux uncertainties derived using our Monte-Carlo simulations are denoted 
by cr s,mu . These flux uncertainties present the advantage of accounting for 
nearly all sources of noise, which explains why they are almost always larger 
than noise estimates based on residual maps (i.e., a Jnap ). However, this noise 
estimate is computed independently of the actual position of the individual 
sources, it is statistical. In some cases, local effects can dominate the noise as it 
is the case when two sources are blended. This local effect, together with the 
background fluctuation due to the photometric confusion noise (i.e. the noise 
due to sources fainter that the detection limit that were not subtracted from 
the image to produce the residual image), is better accounted for in the noise 
estimated from the residual maps, which is estimated locally. To be conservative, 
users should always use the highest uncertainties between a Jnap and cr™'", but 
not the quadratic combination of both since they are not independent. 

Tables B.8, B.9, B.10 and B.ll give excerpt of our complete GOODS-N/S 
24 fim and 70 fim catalogs available at CDS (http://cdsweb.u-strasbg.fr/cgi- 
bin/qcat?I/A+A/<volume>/<page>). For each field we decide to split our 
24 fim catalogs into two (i.e., sources with 3 < cr s """ < 5 and sources with 
5 > cr s,m ") in order to highlight that in deep and confused fields the use of 
sources below 5-cr has to be done with caution. Positions of the 24 fim and 
70 fim sources correspond to the IRAC positions used as priors to our source 
extraction. IRAC coordinates are calibrated to match the GOODS ACS version 
2 coordinate system. For 24 fim sources that are not individually detected at 
70 fim, we report an upper flux limit computed from our residual maps (i.e., 
5-<T map at the position of the source). 
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Table B.l. Parameter values of the rest-frame 8 //m, 15 //m, 25 jim and 35 LF 



Redshift 


Wavelength 


a\ a 


Q-2 a 


Log(L kn „) 


Log(</w) 










Log(L ) 


Log(Mpc~ 3 dex~ l ) 


z~0 


8 /im 


-0.8 


-3.2 


10.29 ± 0.01 


-2.92 ±0.01 


1.3 < z< 1.8 


8 //m 


-0.8 


-3.2 


11.33 ±0.07 


-3.35 ±0.09 


1.8 < z < 2.3 


8 [im 


-0.8 


-3.2 


11.41 ±0.06 


-3.63 ± 0.09 


z~0 


15 //m 


-0.57 


-2.27 


9.56 ± 0.04 


-2.73 ± 0.07 


0.4<z<0.7 b 


15 //m 


-0.57 


-2.27 


10.22 ± 0.03 


-2.63 ± 0.05 


0.7 < z < 1.0 b 


15 ^/m 


-0.57 


-2.27 


1 57 + 04 


-2 86 + 04 


1.0 <z< 1.3 b 


15 //m 


-0.57 


-2.27 


10.79 ± 0.05 


-2.93 ± 0.06 


1.3 < z< 1.8 


15 //m 


-0.57 


-2.27 


10.85 ± 0.06 


-3.02 ± 0.08 


1.8 < z< 2.3 


15 //m 


-0.57 


-2.27 


10.99 ± 0.08 


-3.17 + 0.11 


z~0 


25 //m 


-0.6 


-2.2 


9.91 ±0.01 


-3.04 ± 0.02 


1.3 < z< 1.8 


25 fim 


-0.6 


-2.2 


11.12 + 0.11 


-3.12 + 0.13 


1.8 < z< 2.3 


25 //m 


-0.6 


-2.2 


11.18 ±0.07 


-3.14 ± 0.11 


z~0 


35 //m 


-0.55 


-1.95 


9.85 ± 0.07 


-2.83 ±0.10 


0.4<z<0.7 b 


35 //m 


-0.55 


-1.95 


10.73 ± 0.04 


-2.98 ± 0.05 


0.7 <z< 1.0 b 


35 //m 


-0.55 


-1.95 


10.82 ± 0.04 


-2.73 ± 0.06 


1.0 <z< 1.3" 


35 fim 


-0.55 


-1.95 


11.20 ±0.04 


-2.98 ± 0.05 


1.3 <z< 1.8 


35 //m 


-0.55 


-1.95 


11.26 ± 0.11 


-3.12 ± 0.12 


1.8 < z< 2.3 


35 fim 


-0.55 


-1.95 


11.34 ±0.07 


-3.10 ±0.09 



a Fixed slopes 

b These parameter values are taken from Magnelli et al. (2009) 



Table B.2. The rest-frame 8 /im LF derived from the 1/V max analysis 



1.3 < z< 1 


.8 


1.8 < z< 2.3 




logW 


i°g(^)- lo s(0 


log(0) 


log(ie) 


log(MpcT 3 dex~') 


log(Lo) 


log(Mpc _3 dex _1 ) 


10.4 - 10.6 


-2-68!°:!° 


10.6-10.8 


-3 14 +0U 

mmi 


10.6-10.8 


9 ss+o.io 
_2 - ! "-o.io 


10.8- 11.1 


-3 16 +ou 


10.8- 11.0 


_j Q4+0.10 


11.1 - 11.4 


o 47+0.11 


11.0- 11.2 


t i 9+0.II 
J - LA -0.ll 


11.4- 11.7 


-3 96 +014 


11.2- 11.6 


J -0.12 


11.7 - 12.0 


-5 14+ 031 

J - -5.14 


11.6- 12.0 


-4.63!j5i? 






12.0 - 12.4 


-5.23!°j> 







Table B.3. The rest-frame 15 yum LF derived from the 1/V max analysis 



1.3 <z< 1. 


8 


1.8<z<2.3 


log(L' S ; ra )-log(L'-; m ) 
log(Lo) 


log(0) 
log(Mpc _3 dex _1 ) 


lOg^m)" 1 ^^) 

log(Lo) 


log(0) 
log(Mpc _3 dex _1 ) 


10.1 - 10.3 


_? 54+0.10 
z - J ^-o.n 


10.4- 10.6 


9 77+O.1 1 
' i ' ' '-0.11 


10.3 - 10.6 


9 7Q+0.10 
^•"-o.io 


10.6 - 10.9 


-3 05 +01 ° 


10.6 - 10.9 


-3 or 010 

J - w -o.io 


10.9- 11.2 


-2 07+0.11 
-0.11 


10.9- 11.2 


-3 45 +011 

J -^ J -0.12 


11.2- 11.6 


-4 01 +013 


11.2- 11.7 


_ 4 ^9+0.15 
H - J -0.19 


11.6- 12.1 


-5 06 +0 - 25 

J -" -0.54 


11.7- 12.2 


-5 02 +0 ' 25 
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Fig.B.l. The rest-frame 8 fim LF estimated in two redshift bins with the l/V max method. Red squares are obtained using scaled CE01 templates 
which best fit the L^m/ci+z) _ i-70/jm/a+z) correlation. Empty triangles and blue dashed-dotted line present the rest-frame 8 jim LF obtained at 
z ~ 2 by Caputi et al. (2007). Asterisks show the local reference taken from Huang et al. (2007) and the dotted line presents the best-fit to these 
data points with a double power law function with fixed slopes (see Table B.l). The dark shaded area span all the solutions obtained with fhe;f 2 
minimization method and compatible, within 1 cr, with our data points. The dashed line represents the best fit of the rest-frame 8 jim. LF. 
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Fig.B.2. The rest-frame 15 yum LF estimated in two redshift bins with the \IV max method. Red squares are obtained using scaled CE01 templates 
which best fit the L24^m/(i+z) _ ^70^m/(i+z) correlation. Asterisks show the local reference taken from Xu (2000) and the dotted line presents the 
best-fit to these data points with a double power law function with fixed slopes (see Table B.l). The dark shaded area span all the solutions obtained 
with the^ 2 minimization method and compatible, within 1 cr, with our data points. The dashed line represents the best fit of the rest-frame 15 
yum LF. In the first redshift panel, we reproduce in green, blue, yellow and red the best fit of the LF obtained at 0.4 < z < 0.7, 0.7 < z < 1.0, 
1.0 < z < 1.3 (Magnelli et al. 2009), and 1.8 < z < 2.3 respectively. 
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Fig. B.3. The rest-frame 25 yum LF estimated in two redshift bins with the 1/V„, av method. Red squares are obtained using scaled CE01 templates 
which best fit the L^m/ri+z) - Ljo^m/o+z) correlation. Asterisks show the local reference taken from Shupe et al. (1998) and the dotted line presents 
the best-fit to these data points with a double power law function with fixed slopes (see Table B.l). The dark shaded area span all the solutions 
obtained with the^ 2 minimization method and compatible, within 1 cr, with our data points. The dashed line represents the best fit of the rest-frame 
25 fim LF. 




Fig. B.4. The rest-frame 35 yum LF estimated in two redshift bins with the 1/V„,„ V method. Red squares are obtained using scaled CE01 templates 
which best fit the i.24/jm/(i+z) ~ ^70;im/(i+z) correlation. Asterisks show the local reference derived from Shupe et al. (1998) and the dotted line 
presents the best-fit to these data points with a double power law function with fixed slopes (see Table B.l). The dark shaded area span all the 
solutions obtained with the x 1 minimization method and compatible, within 1 cr, with our data points. The dashed line represents the best fit of 
the rest-frame 35 fim LF. In the first redshift panel, we reproduce in green, blue, yellow and red the best fit of the LF obtained at 0.4 < z < 0.7, 
0.7 < z < 1.0, 1.0 < z < 1.3 (Magnelli et al. 2009), and 1.8 < z < 2.3 respectively. 
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Table B.4. The rest-frame 25 fim LF derived from the 1/V max analysis 



1.3 < z < 1. 


8 


1.8<z<2.3 


loe(L'°' v ) - los.(L"' s " ) 

lu t,y^25iim' lu &^25/im' 1 
log(Zo) 


log(Mpcr 3 dex _1 ) 


iog(4 o 5 ; m )-iog(^ m ) 

log(Lo) 


log(</>) 
log(Mpc _3 dex _1 ) 


10.3 - 10.6 


— / n5 

z - OJ -o.io 


10.6- 11.0 


-2 90 +01 ° 

■""-o.io 


10.6- 10.9 


9 09+O.IO 

~ 2 - s/ -o.io 


11.0- 11.4 


q 17+0.10 
J,i -0.10 


10.9- 11.2 


t 91+O.II 

J - -o.n 


11.4- 11.8 


-4 06+ 013 

^•""-0.16 


11.2- 11.7 


74+0.11 
J - ' -0.12 


11.8- 12.2 


-4 Q7+0.25 


11.7- 12.2 


-4 85 +0 - 22 

° -0.38 







Table B.5. The rest-frame 35 /im LF derived from the 1/V max analysis 



1.3 <z< 1. 


8 


1.8<z<2.3 


log(^)-log(%) 

log(Zo) 


log(Mpc _3 dex _1 ) 


log(4- m )-log(L^) 
log(Lo) 


log(0) 
log(Mpc _3 dex _1 ) 


10.5 - 10.8 


9 7O+0.10 

" 2 - /8 -o.io 


10.7 - 11.0 


-2-8011! 


10.8 - 11.1 


9 qo+0.10 


11.0- 11.2 


9 CH+0.11 
^-^-o.n 


11.1 - 11.4 


-3.22^1 ! 


11.2- 11.6 


-3 16 +0 '° 

J - 1 -0.10 


11.4- 11.9 


-3.76!»;! 2 


11.6- 12.1 


-3 96+ 012 

-0.13 


11.9- 12.4 


-4.85!»;| 


12.1 - 12.6 


-5 36 +031 

J - J -5.36 



Table B.6. The infrared LF derived from the l/V max analysis using L 



fit 



1.3 < z< l.S 



log(L;-)-log(L7 ! ) 
log(Lo) 



1.8 < z < 2.3 

log(Lf/) - log^f) 



log(0) 
log(Mpc _3 dex _1 ) 



log(Lo) 



log(0 
log(Mpc _3 dex _1 ) 



10.8- 


11.2 


-2-38l : | 


11.25 - 11.55 


9 79+O.25 
' -0.25 


11.2- 


11.6 


-2.78H 


11.55 - 11.85 


-3 05 +0 25 

-0.25 


11.6- 


12.0 


-3.1511 


11.85 - 12.15 




12.0- 


12.4 


-3 69 +0 - 26 


12.15 - 12.45 


1 oq+0.26 


12.4- 


12.8 


-4 75+ - 31 

H -' J -0.45 


12.45- 12.75 
12.75 - 13.05 


-4 84 +0 - 34 

H -° -0.58 

-5 14 +039 



Table B.7. The infrared LF derived from the 1/V max analysis using L]\ 



1.3 <z< 1. 



1.8<z<2.3 



log(Lf/) - log^f ) 
log(Lo) 



log(L^)-log(L;7 ! ) 
log(Lo) 



log(0) 
log(Mpc _3 dex _1 ) 



log(0) 
log(Mpc _3 dex _1 ) 



10.8 - 11.4 
11.4- 12.0 
12.0- 12.6 
12.6 - 13.2 



9 ce+0.25 

-3.1111 
-3-9111 
-5.40l 3 



11.15 - 11.45 
11.45 - 11.75 
11.75 - 12.05 
12.05 - 12.35 
12.35 - 12.65 



_9 74+0.25 
'^-0.25 

_t in+0.25 
J - -0.25 

_■} i 1+0.25 
J.ll_ 025 

-3 76 +a2S 

- 3 - /o -0.26 

-4 84+034 

^•°^-0.58 
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Table B.8. MIPS sources in GOODS-N with S 24pm/o"l'™ m > 5. The second and third columns list the prior position from IRAC. Fourth column 
lists the 24 /im flux density. Fifth and sixth columns give the 24 /im flux uncertainty derived from the residual map and from our Monte-Carlo 
simulations respectively (see text for detail). Seventh column indicates the integration time on the source. Columns 8 to 1 1 repeat columns 4-7 for 
70 urn. 



Name 




IRAC position 




Map 
u 24pm 


u 24pm 






Map 
u 70pm 


rr Simu 
" 70pm 


Covn )llm 






anooo i5j2ooo 


lily 


//Jy 


//Jy 


s 


mJy 


mJy 


mJy 


s 


MIPSJ123539.5+621 129.0 


12 


35:39.48 +62:11:29.02 


42.6 


5.6 


6.7 


10873 


< 2.9 






5601 


MIPSJ123539.5+621243.8 


12 


35:39.54 +62:12:43.78 


54.0 


4.6 


6.7 


10239 


< 2.2 






4747 


MIPSJ123539.9+621324.8 


12 


35:39.94 +62:13:24.79 


52.2 


4.2 


6.7 


9512 


< 2.4 






3728 


MIPSJ123540.2+621224.2 


12 


35:40.17 +62:12:24.16 


107.1 


7.2 


8.2 


11518 


< 3.5 






5453 


MIPSJ123540.7+621218.9 


12 


35:40.70 +62:12:18.93 


107.1 


5.8 


8.2 


11934 


< 3.5 






5711 


MIPSJ123541.0+621136.1 


12 


35:41.00 +62:11:36.13 


139.8 


4.1 


8.3 


12437 


< 4.6 






6409 


MIPSJ123541.4+621217.4 


12 


35:41.39 +62:12:17.37 


419.5 


5.2 


11.7 


12395 


2.6 


0.3 


0.6 


6102 


MIPSJ123541.4+621316.5 


12 


35:41.41 +62:13:16.51 


41.3 


4.1 


6.6 


12733 


< 1.6 
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MIPSJ123541.6+621151.0 


12 


35:41.60 +62:11:51.01 


44.7 


3.5 


6.7 
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<4.2 
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MIPSJ123541.7+621223.5 


12 


35:41.73 +62:12:23.54 


64.0 


5.2 


7.0 


12853 


< 2.4 
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Table B.9. MIPS sources in GOODS-N with 3 < S 2 4vJ°~fZ m < 5 - Columns are the same as in Table B.8. 



Name 




IRAC position 


FlS,p. m 


Map 
u 24pm 


_S imu 
u 24pm 


C()V2Apm 


FlOpm 


Map 
" 70pm 


rr Simu 
u 70pm 


CoVjOfim 
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lily 


lily 


s 
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mJy 


mJy 


S 


MIPSI123540.2+621 108.2 


12 


35:40.23 +62:11:08.16 


23.8 


4.2 
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< 2.0 
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12 


35:41.29 +62:10:47.22 
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3.7 
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< 3.4 
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MIPSJ123543. 8+621218.8 


12 
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23.5 


3.5 


6.8 
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< 2.4 






7117 


MIPSJ123544.7+621246.8 


12 


35:44.69 +62:12:46.78 


20.5 


5.7 


6.7 
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< 1.3 






6439 


MIPSJ123545.2+621 151.9 


12 


35:45.18 +62:11:51.91 


23.3 


5.4 


6.8 


15912 


< 3.0 






7764 


MIPSJ123545. 3+621 134.0 


12 


35:45.30 +62:11:33.97 


20.5 


4.5 


6.7 


16640 


< 4.6 






7482 


MIPSJ123545.4+621306.4 


12 


35:45.44 +62:13:06.39 


31.6 


3.7 


6.6 


18447 


< 1.5 






6228 


MIPSJ123545.6+621034.7 


12 


35:45.63 +62:10:34.68 


29.3 


7.1 


6.6 


13190 


< 3.1 






6388 


MIPSJ123547.6+621 147.2 


12 


35:47.59 +62:11:47.18 


31.2 


7.2 


6.7 


18975 


< 3.5 






8395 


MIPSJ123550.3+621423.6 


12 


35:50.33 +62:14:23.64 


21.7 


4.2 


6.8 


13135 


<3.8 






4727 



Table B.10. MIPS sources in GOODS-S with S^pm/v*™" > 5. Columns are the same as in Table B.8. 



Name 






IRAC position 


F2Apm 


Map 
u 24pm 


u 24pm 


C()V 2 4pm 


FlOpm 


Map 
u 70pm 
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<Sj2000 
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S 
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s 
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1 


274331.1 


03 


32:01.11 


-27:43:31.07 
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3.7 
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17558 
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7946 


MIPSI033201 


2 


274636.0 


03 


32:01.15 


-27:46:35.98 


200.7 


6.2 


8.1 


10241 


< 1.7 






8105 


MIPSI033201 


2 


274134.6 


03 


32:01.24 


-27:41:34.57 


186.0 


5.8 


7.7 


12533 


< 2.3 






6676 


MIPSI033201 


3 


274553.8 


03 


32:01.26 


-27:45:53.75 


72.3 


4.4 


7.4 


12171 


< 2.2 






8099 


MIPSJ033201 


4 


274646.5 


03 


32:01.40 


-27:46:46.54 


167.3 


8.0 


7.9 


10220 


< 1.4 
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5 


274138.7 


03 


32:01.46 


-27:41:38.71 


208.1 


6.0 


8.2 
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2.8 


0.2 


0.8 


6689 
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5 


274229.8 


03 


32:01.48 


-27:42:29.80 


95.7 


4.5 


8.1 


20170 


< 5.5 






7206 


MIPSJ033201 


5 


274402.5 


03 


32:01.51 


-27:44:02.45 


57.1 


4.3 


6.7 


17091 


< 2.3 






8299 


MIPSJ033201 


6 


274326.9 


03 


32:01.61 


-27:43:26.93 


55.4 


2.8 


6.8 


19023 


< 1.8 






7935 


MIPSJ033201 


7 


274349.4 


03 


32:01.66 


-27:43:49.37 


51.0 


4.1 


6.7 


18518 


< 1.8 
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Table B.ll. MIPS sources in GOODS-S with 3 < S 2 4pm/cr^ m < 5. Columns are the same as in Table B.8. 



Name 


IRAC position 


F24pm 


Map 
(T 24pm 
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CoV24pm 
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Map 
°~70pm 


rr S imu 
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5J2000 
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s 


MIPSI033158.1-274207.9 


03:31:58.08 


-27:42:07.92 


23.5 


5.0 


6.8 
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< 2.3 
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03:31:58.93 


-27:43:59.09 


27.2 


5.3 


6.8 


10417 


< 1.5 






7976 


MIPSJ033159.1-274421.1 


03:31:59.12 


-27:44:21.14 


29.8 


4.6 


6.7 


9963 


< 2.1 






8036 


MIPSJ033159.2-274145.8 


03:31:59.18 


-27:41:45.82 


22.5 


6.9 


6.9 


10825 


< 1.5 






6791 
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03:31:59.71 


-27:44:59.82 


21.7 


6.8 
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< 2.9 






8135 
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-27:45:26.57 


20.9 
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5.6 


6.6 
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27.0 
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22.9 


4.6 
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15533 
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